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ABSTRACT 

Supermassive black hole (BH) masses (A^bh) are strongly correlated with galaxy stellar bulge masses 
(A^buigc) and there are several ideas to explain the origin of this relationship. This study isolates 
the role of galaxy mergers from considerations of other detailed physics to more clearly show how a 
linear ^Yi-galaxy mass relation (A^BH--^gai) can naturally emerge regardless of how primordial BHs 
were seeded inside galaxies, if the galaxy mass function declines with increasing mass. Under this 
circumstance, the A^BH-A^gai relation is a passive attractor that eventually converges to a tight linear 
relation because of two basic statistical effects: a central limit-like tendency for galaxy mergers which 
is much stronger for major mergers than for minor mergers, and a convergence towards a linear relation 
that is due mainly to minor mergers. A curious consequence of this thought experiment is that, if 
galaxy bulges are formed by major mergers, then merging statistics naturally show that A^bh would 
correlate more strongly with bulge dominated galaxies, because of stronger central-seeking tendencies, 
than with disk dominated galaxies. Even if some other physics is ultimately responsible for causing 
a linear AlBH-Albuigc relationship, this thought experiment shows that, counter to intuition, random 
merging of galaxies tends to strengthen rather than weaken a pre-existing, linear, correlation. This 
idea may be generalized to other gravitationally bound systems (dark matter halo, compact nuclear 
objects) that retain their physical identities after experiencing mergers. 

Subject headings: galaxies: bulges — galaxies: formation — galaxies: evolution - galaxies: statistics 



1. INTRODUCTION 

In recent years, there have been several surpris- 
ing discoveries of fundamental scaling relations be- 
tween supermassive black hole masses (A^bh) and 
large scale galaxy bulge properties: stellar veloc- 
ity dispersion cr*, bulge mass bulge, the pro- 
file slope of galax i es, a nd the inner core ra- 
dius (iGebhardt e t all 120001: iFerrarese fc Merrita 120001: 



Kormcndv & Rich stond 119951 iMagorrian et al.lll998l: iHol 
^99; Kormcndv fc Gebhardtl^2001': 'Graham et al .ll200ir 



Marconi & Hunt 2003; Haring & Rix 2004; Bart h et al 
2005 t iLauer et al.1 12006t iGreene fc Hoi 120061: IWoo et al 
20061 ). At masses lower than A^buige^ 10^°A^n, the cen- 



tral object with which galaxies correlat e may be either in- 
term e diate mass BHs in dwarf g a laxies ([Filippeiiko fc Hoi 
[200l iBarth et all l200l . [20051 : IGreene fc Hoi 120061) or 
central massive star clusters ( Ferrarese et al .1 120061 : 
I Graham fc Driverl 120071 ) . The small amount of intrin- 
sic scatter between black holes and bulges is often inter- 
preted to suggest a causal connection between the two 
— that th e growth of one inight somehow regulate the 
other (e.g. lSilk fc R ees' 19981: IPi Matteo et ahllloOS '). 

When did the fundamental BH scaling relationships 
come about? At higher redshifts, observations are still 
in the early stages, complicated in part by the chal- 
lenges of measuring the BH mass and the bulge prop- 
erties in the same galaxy. However, recent evidence 
from the study of quasar host galaxies indicates that a 
fundamental A^BH-A^buigo correlation might h ave been 
present as early as z ~ 4 ijPeng et all l2006allbl ). Fur- 
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thermore, there also appears to be an evolution in the 
AlBH-Albuige ratio by a factor of ~ 4 in the same studies, 
which points to the possibility that the BH masses may 
have matured more quickly than their surrounding stellar 
bulge mass in the past. In addition, other observations 
that use [Oili] and CO emission line widt hs to infer the 
gravi t atiori a ,! potential of quasar bulges (IShields et al.l 
l2(M [20061 : iSalviander et al.l[2(M 120071 : lHoll2007f ) sug- 
gest similar trends^, and residual traces o f evolution re- 
main detectable even below z = 1 (jTreu et a l. 2004, 2007|; 
I Woo et all l2006l ). De spite the general agreements ob - 
servationally (but, see lLi et al.ll2006l : iBorvs et al.ll2005l ). 
there remain several weaknesses that still complicate the 
interpretation: a significant factor of 2 systematic un- 
certainty in the normalization of the BH mass calibra- 
tion, the possibilit y that the normalization of the virial 
BH mass es timate ("Kaspi et al.ll2000!; 'Vestergaard"2002| 



Onken et al. 2004: Kaspi et al. 2005; Greene fc Ho 200i 



Peterson fc Bentzl 120061 : IVestcrgaar d fc Petersoni l200l l 



may evolve with time, because of the fact that the bulge 
masses are not directly measured, and the possibility 
that the host g alaxy mass may b e biased low in high 
redshift quasars (jLauer et al.ll2007l ) because of the steep 
decline in the luminosity function of galaxies. With re- 
gard to the latter selection bias, it is worthwhile to note, 
ho wever, that low redshift {z < 0.3) quasars host galax- 
ies (jMcLure fc Dunlop||200lD . which have si milarly high 
■ Mbh a nd A^gai to the high redshift objects in lPeng et al.1 
(|2006al|bl ). do not exhibit a mass dependent bias even 
when the scatter is large. Despite these uncertainties, 
the finding of an existing strong A^BH-A^buigc correla- 

A lack of evolution in the stellar velocity dispersion im- 
plies that bulge mass decreases with look-back time (see, e.g, 
IRobertson et al.|[2006bl V 
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tion at 2 > 1 is likely to be more secure. 

Controversies about the evolution aside, it is not only 
puzzling that the A^BH-bulge correlations should ex- 
ist, but that they would have a small intrinsic scatter 
of 0.3 dex in A^bh and that t he correl ation wit h the 
bulge mass is pract i cally linear (iMarconi fc Hunt l2003l : 
iHaring fc Rb3 120041 : iLauer et all I2006D . These curious 
facts have received wide theoretical attention over the 
years and could be explained in a number of ways. For 
instance, both the A^BH-A^buigc and AIbh-o'* relation 
can be pro duced by g as accretio n onto a nuclear disk 
(jBurkert fc Silk 2001; Cenl l2007t l or turbulent dissipa- 
tion of gas ((Escala, 200(1 ) followed by star formation, by 
a combinat ion of BH accretion, galaxy mergers, and star 
formation (Li et al. 2006'; Kauffmann fc Haehneltl l2000l : 
iHaehnelt^ Kauffmann 2000), by gravitational collapse 
of inner parts of a galaxy t hat forms the bulge from 
a rotating isother mal sphere (lAdams et al.ll2001l I2003D . 
by stellar capture fMir alda-Escude fc Kollmeieill2005l ) in 
the accretion disk. The effect of dissipation in galaxy 
merging on the BH and fun damental plane scaling rela - 
tions has also b een visited: ICiotti fc van Albadal ()200lD 
speculated, and lRobertson et al.l I 2006af ) show, that dis- 
sipational mergers can produce the A^bh-c* slope and 
maintain the fundamental plane of elliptical galaxies. On 
the other hand, whereas dissipation/ess mergers of ellipti- 
cal galaxie s also tend to preserve the ob s erved fundamen- 
tal pla ne (|Bovlan-Kolchin et al.l 120061 : [Robertson et al.l 
l2006al ). dissipationless mergers of disk galaxies tend to 
produce a relation that more closely parallels the virial 
plane. 

Despite the aforementioned models, the explanation 
that has spawned immense activity is the theory of feed- 
back from an active galactic nucleus (AGN). The AGN 
feedback idea rests observations that quasars typically 
radiate at a fixed fraction (10%— 100%) of the Edding- 
ton luminosity. If such radiation is produced by a massive 
enough BH, this energy budget is in principle sufficient 
to quench s tar formation and t e rminate the BH growth 
itself (e.g. 'Silk fc Rees' Togs'; 'Di Matteo et al.' '2005'; 
Cattanco et al. 2005; Springcl ct al. 2005; Hopkins et al. 
20071 12006al . and references therein) . This idea has been 
incorporated into cosmological merger simulations (e.g. 



Granato et al.ll2004 iFontanot et"alTl2006t ICroton et aD 
20061), and is seen to have profound possibilities for ex- 
plaining a wide array of other galaxy evolution puzzles, 
including the evolution in the galaxy and quasar luminos- 
ity functions, mass functions, star formation rates, the 
X-ray b ackground, and the b imodality of galaxy colors 
(e.g. see 'Hopk ins et aLll2006al lbl. and references therein). 

Even though AGN feedback is promising for explaining 
many aspects of galaxy evolution, it is possible that some 
other mechanism (or mechanisms) may either have signif- 
icant influence on the scaling relation between A^bh and 
A(buigcj or in fact be more fundamentally the cause. It is 
therefore worthwhile to look for such factors and to fully 
study their effects. Dating back years before AGN feed- 
back physics became a popular idea, one such fundamen- 
tal factor considered by many was galaxy merging (e.g. 
Kauffmann fc Haehneltl '200 0|: | Haehnelt fc K auffmanri 
2D00: Ciotti fc van Albada~ l20pil: iNipoti e~ al. 2003; 
Islam et al.l I2003L 120041 : IVolonteri et all I2003IV Even 



fering degrees of agreement with the observations. How- 
ever, the reason why a tight linear correlation should 
emerge is not immediately apparent even when consider- 
ing no other physics besides merging. In fact, as accorded 
correctly by intuition, a correlation indeed cannot arise 
spontaneously from chaotic combinations of galaxies and 
black hole masses in general. As will be shown below, 
what does make a tight correlation emerge is the fact 
that the galaxy mass function decreases with increasing 
mass. This circumstance brings about a number of inter- 
esting implications that will be explored in future stud- 
ies. In this study, the modest goal is to show that when 
the focus is switched from the BH-bulge relation, to un- 
derstanding the more general ^Yi-galaxy relation, some 
insights may be gained into understanding the growth 
and evolution of the TWbh-AI bulge relation itself. 

This study therefore revisits the issue of galaxy merg- 
ing from the standpoint of a thought experiment. This 
toy model identifies three basic reasons for why the 
A^BH-A^buigc relation may appear the way it does, even 
if the BH masses and their host galaxy masses were 
completely uncorrelated initially, or if they started out 
with a steep powerlaw correlation. The ov er-arching 
premi se is that the galaxy mass function has a lSchechtei 
powerlaw form, so that there is a decline in galaxy 
number density with increasing mass, especially above 
masses M.* . When this circumstance is met, it can 
be shown that a "linear attractor" and a central limit- 
like tendency can work efficiently to produce a linear 
A^BH-Atgai relation, and to reduce the scatter over time. 
However, while it is tempting to generalize this result 
to the A^BH-f* relation, it is not as simple to do be- 
cause it is not clear how the stellar velocity dispersion 
scales during galaxy mergers, something that depends on 
physics not considered in this study (e.g. dissipational 
vs. dissipationless mergers, star formation, AGN feed- 
back - see Robertson et al.l 12006 b: Hopki ns etlil] 120071 : 
iDi Matteo et al..,2007; .Siiacki et al... 20071 

The main emphasis of the current study is thus 
only to present a pure statistical exercise, and as such 
will not invoke merger trees or external physics - be- 
sides which many and much mo re sophisticated mod- 
els have aheady been run (e.g. Volonte ri et al.l [20031; 
Islam et all 120041: iGranato et al.l l2004 iF ontanot et all 



20061: [Robertson et al.l l2006bl : ICiotti et al.ii2007i) . In a 



without feedback, those simulations already find it pos- 
sible to produce the BH scaling relations, albeit to dif- 



sense, this work examines a common thread shared by 
all previous cosmological simulations. While it is tempt- 
ing to invoke realism by introducing detailed physics 
from the get-go, e.g. star formation, BH accretion, and 
AGN feedback, isolating the effects of simple statistics 
enables a cleaner exposition of why the convergence be- 
havior should be expected. As such, the current study 
is not a critique of other, more physical, models which 
can also explain the same correlation, or to pass criti- 
cal judgment about which is more or less relevant. In- 
stead, the main message of this study is that, regard- 
less of what other physics may ultimately produce the 
AleH-Albuigc relation or weaken it, an existing corre- 
lation should strengthen if galaxies continue to merge 
thereafter, whether by major or by minor mergers. In 
other words, merging alone, in the absence of all other 
physics, is a sufficient condition to bring about a tight, 
linear, A^BH-A^buige relation over time, and is always 
"pulling behind the scenes" to bring about such a corre- 
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lation. The issue of interest to follow up is to what extent 
this, or other as yet identified mechanisms, may matter 
in the end, and to predict more realistic scatter in the 
A^BH-Algai relation (as opposed to only AleH-Albuigc) 
under this hypothesis using more realistic merger trees 
and physics. We will address this subject in a foUowup 
study. 

The following discussion will begin by presenting a 
heuristic view to explain why a linear A^BH-A^gai re- 
lation is a natural consequence of random merging (Sec- 
tion 2), followed by Monte-Carlo simulations to illustrate 
the effect (Section 3), and lastly by a discussion and con- 
clusion. In much of the discussion below, the relationship 
under consideration is more generally the AlBH-A^gai re- 
lation, whereas BHs are thought to correlate more tightly 
with the spheroid component of A^gai, that is bulge- It 
will be seen that the tighter correlation between Mbb 
and A^buige is a special case of the A^BH--A^gai relation 
and can be understood in the same framework if this hy- 
pothesis represents the dominant mechanism by which 
BHs correlate with galaxy masses. 

2. HEURISTIC PICTURES 

2.1. The Central Limit Theorem of Galaxy Mergers 

As galaxies undergo merging, it can be shown that 
the scatter in the A^BH-Algai relation diminishes with 
increasing number of mergers as a consequence of the 
central limit theorem (see Appendix). To see this most 
easily, consider first major mergers, which by definition 
occur between galaxies of roughly equal mass, often de- 
fined to be within the range M1/M2 < 4 (Figure [1^). 
Their similarity in mass also means that the number dis- 
tribution (i.e. mass function) of A^bh at a specific galaxy 
mass, Ml, is similar to that at A^2, i-e. they are drawn 
from parent distributions that have roughly the same 
shape (Figure [1)3) . However, the mean (/i) of the distri- 
butions might be offset by an amount that depends on 
the steepness of the A^BH-Algai correlation: the steeper 
the correlation (Figure [T^, upper ellipse) , the larger is 
the offset. 

Furthermore, consider that in the limit of no initial 
A^BH--A^gai correlation (Figure lower ellipse), /i is 
constant with galaxy mass. The average of any two 
BH masses randomly drawn from the parent distribu- 
tion (e.g. Fig. [T|d) therefore has a tendency toward 
the mean value of the parent, by the central limit 



theorem. This tendency also means that the resulting 
BH mass distribution from mergers, obtained by sum- 
ming BHs drawn from the same distribution, will have a 
smaller log-normal scatter, cr(log(/i)) = cr(/i)//i, than the 
original log-normal distribution, because the fractional 
mass error, a {Mbk,i + Mbk,2) / (A^BH.i + A^bh,2), 
is smaller, which means the scatter in the A^bh vs. 
Algal relation, logarithmic on both axes, decreases. 
In the instance where the initial AIbh mass distribu- 
tion (Figure [1)3) is normal, with a scatter (TBH.init, the 
resulting scatter of all galaxies that have undergone 
one full cycle of major mergers is a (log /j.bh, merge) — 
cr (log /UBH.init) /V2 (see Appendix). Therefore an en- 
semble of galaxies that has undergone A'maj mergers 
should have a scatter in the AJbh- Algal relation that 
is reduced by ~ 2°'^^"='-i , compared with the initial rela- 
tion. However, a log-normal AIbh distribution will have 
a different convergence rate. 

Clearly, this central limit theorem behavior applies to 
a finite parent distribution of any shape, but the size of 
the scatter and the rate at which the scatter decreases 
both depend on the shape of the distribution and the 
steepness of the initial AIbh- Algal correlation. In the 
situation where the correlation between AIbh and A^gai 
is steep (e.g. Fig. [T^, upper eUipse) the effective cu- 
mulative mass function of BHs (Figure lb) residing in 
galaxies with J\42 ^ Algal ^ All tias a wider (TBH.init 
than if the AlBH-Algai relation is shallow (Fig.[T^, lower 
ellipse). For the same reason, minor mergers also have 
wider (TBH.init, as the galaxy mass differences are larger, 
than do major mergers. Therefore, galaxies that have 
only undergone major mergers would produce an AIbh- 
Atgai relation that is significantly tighter than for galax- 
ies that have only experienced minor mergers, for the 
same merger rate, and if the initial mass correlation is 
not flat. This effect is seen in the Monte-Carlo simula- 
tions below. 

2.2. Galaxy Merging From a Schechter Mass 
Distribution 

If the mas s density of gala xies follows a Schechter pow- 
erlaw form (|Schechteilll976D . 

^(^) = 'J>o(;^)""oxp(--^), (1) 
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log (Mgai ) log (Mgai ) 



Fig. 2. — (a) No initial correlation in the jMeH-'^gal relation. A galaxy at location 1 is more likely to merge with another galaxy 
at a much lower A^gai than with one that is comparable to itself (due to the bottom heavy mass function), but roughly of comparable 
-Mbh (due to non-correlation between jVfBH-A^gal)- Thus the net evolutionary vector is steep. In contrast, a galaxy at location 2 is likely 
to merge with another that is comparable in both Mbh ^^nd Algal to itself, so the evolution vector is shallower. (6) An initial, strong 
powcrlaw, correlation in the A^BH-A4gai relation. A galaxy at location 1 is more likely to merge with another galaxy comparable in A^gai 
but at a much lower A^bh than itself. Thus, the net evolutionary vector is shallow. In contrast, a galaxy at location 2 is likely to merge 
with another of comparable in AIbh and Algal to itself, so the evolution vector is steeper. 



then it is possible also to show that a linear A^BH--A^gai 
relation naturally emerges over time, so that the relation, 

XbhW =r(z)Mgai(z)'3, (2) 

eventually takes on /3 = 1. The value of F, which lo- 
cally is measured t o be r(0) ~ 1/800 for bulges (e.g. 
iHaring fc Rixll20(j^ . is otherwise arbitrary in the discus- 
sion below, r is degenerate with respect to assumptions 
about the initial scatter of the A^BH-A^gai relation, the 
initial slope /3, and the initial normalization r(oo), for 
which there are currently insufficient observational con- 
straints; it will not be addressed further in this study. 
The other assumption used here is that the probability 
for two galaxies to merge comes from Monte-Carlo sam- 
pling of a Schechter mass function (Eq. 1). In actuality, 
galaxies do not merge randomly, especially at late times. 
However, complete randomness is only used to facilitate 
the discussion, and is not a pre-requisite, since the rea- 
soning depends only on the fact that minor mergers oc- 
cur more frequently than major mergers. This assump- 
tion does mean that, depending on the relative balance 
of major vs. minor mergers, the effects described here, 
namely convergence toward linearity versus central-limit 
behavior, may be more relevant at some epoch in time 
than at others 

The reason that a linear correlation emerges through 
galaxy mergers is illustrated in Figure 2. Figure 2a shows 
the situation in which the initial BH and galaxy mass 
distributions are completely uncorrclated, so that (3 = 0. 
The lower part of the diagram shows hypothetical mass 
functions with two different "faint end" slopes, A and B, 
which will be individually considered in the Monte-Carlo 



simulation below. If there is no correlation between the 
A^BH and A^gai, then the ratio A^Bn/Algai will be, on 
average, larger for low mass galaxies than for high mass 
galaxies as can be seen by comparing the A^Bu/A^gai ra- 
tio at any two locations, for example, those labeled "1" 
and "2." Therefore, as galaxies merge, a massive object 
at the extreme end of the mass function, located at posi- 
tion 1, 071 average, is more likely to merge with another 
having a much larger A^BH/-A^gai, thereby evolving the 
merger product in a steep upward direction, as exagger- 
ated by the vertical arrow. In contrast, an object at 
position 2 is likely to merge with objects comparable in 
both A^BH and A^gai to itself, so the net evolutionary 
vector has a shallower slope. Therefore, the cumulative 
effect of mergers along the mass spectrum is to steepen 
the massive end of the A^BH-A^gai relation relative to 
the lower extreme, even as the lower end grows in A^bh 
and A^gai on average. 

In the other extreme. Figure 2b, if the primordial rela- 
tion between the BH and galaxy masses is steep, corre- 
sponding to (3 > 1, the opposite behavior occurs. Galax- 
ies at location 1 generally have a larger A^BH/A^gai ra- 
tio than galaxies that have lower mass. Therefore, the 
A^BH/A^gai ratio for massive galaxies would tend to de- 
crease through mergers. The net effect on massive galax- 
ies is to evolve the merger remnant more quickly to the 
right on average than a lower mass galaxy at location 2. 

Because of a mirror symmetry between Figures 2a and 
2b, the natural equihbrium state of the A^BH-A4gai re- 
lation is at /3 = 1, so that further merging of galaxies 
would evolve remnants along the linear relation, with a 
constant ratio r=A^BH/A/'gai- Also because of the con- 
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vergence toward this "attractor" state the scatter in the 
relationship would necessarily decrease over time through 
galaxy merging. 

Lastly, it is worth mentioning that the presence of a 
break in the galaxy mass function at M* is a sufficient, 
but not necessary, condition for convergence toward lin- 
earity. A pure powerlaw with A4*— co would produce a 
similar behavior, however, the convergence is slower for 
flatter slopes (a — > —1). Furthermore, while the con- 
vergence behaviors just described is quite strong for a 
Schechter mass function with a = —I, the convergence 
would fail for a pure powerlaw with the same slope, be- 
cause of a lack of a break in the mass function. 

3. MONTE CARLO SIMULATIONS 

To illustrate the idea discussed above, and to quantify 
how quickly a linear A^BH-A^gai relation might emerge, 
it is useful to consider several numerical simulations for 
the situations shown in Figure 2. For each of the two 
scenarios, no-correlation (Fig. 2a) and steep correlation 
(Fig. 2b), it is also instructive to consider two differ- 
ent initial mass functions, A and B, shown in the lower 
half of Figure 2. The two powerlaw slopes explored be- 
low are a = —1.5 and a = —0.5, respectively. These 
choices are motivated by observations of the lum inosity 
funct i ons for high rcd shift galaxies (e.g. Gabasch^et_alJ 
120061 : iGiallongo et al.l l2005h under the assumption that 
light traces mass. 

3.1. Simulation Set-up and Definitions 

Definition of the number of major and minor mergers. 
— As implied in Section 2, how closely a galaxy lies 
to the linear part of the A^BH-A^gai relation, and how 
tight the final scatter is for an ensemble of galaxies, will 
depend on the cumulative merger history. Therefore, the 
most useful way to understand the simulation results is 
to define the number of major mergers, A'maj, as the 
cumulative sum of all such events over the entire tree for 
a given galaxy, not just in the most immediate, that is, 
the main, branch. For example, even if a galaxy has never 
experienced a single major merger on the main branch 
in its lifetime, it could still lie close to the final, linear 
part of the A^BH--A^gai relation because the progenitors 
of the main branch, and their progenitors, and so forth, 
could have experienced a number of major mergers. 

Evolution of the mass function. — One issue to con- 
sider is how the galaxy mass function might evolve, and 
whether the path of evolution might affect the final con- 
vergence. The effect of galaxy masses growing with time 
is to both increase Ai* and steepen the "faint end slope" 
(a) of the Schechter function. As galaxy populations 
grow in mass and number density, the rate of change in 
A4* and a would affect the rate of convergence to the 
final AlBH-A^gai normalization, slope, and scatter. Pre- 
dicting the rate of change in the A^BH-A^gai relation re- 
quires realistic merger trees and accounting for other de- 
tailed physics such as feedback, which will be addressed 
in a foUowup study. For the current purposes of show- 
ing that convergence toward a linear AlBH--A^gai relation 
does naturally occur, it suffices to consider two scenarios: 
a replenishment scenario, in which the Schechter mass 
function is continuously, and randomly, replenished as 
galaxies merge, and a depletion scenario, in which no new 
galaxies are formed to take the place of those that have 



merged. Combined with considerations about the initial 
mass function slopes, the simulations will have covered 
the gamut of sensible possibilities and conditions that 
may be present at early and late cosmic epochs. 

Initial scatter of the A^BH-A^gai relation and the initial 
mass function M* . — In the simulation, the distribution 
of galaxies is first drawn randomly from the Schechter 
mass function, after which a BH mass is assigned, fol- 
lowing Equation 2, by drawing from a log- normal Gaus- 
sian distribution with a generous Gaussian dispersion of 
(Jbh — 2 dex, i.e. a scatter of a factor 100 in mass. 
The exact choice of the initial scatter is directly propor- 
tional to, but otherwise only partially determines, the 
final scatter in the A^BH--A^gai relation. Other factors 
that determine the final scatter depend on how long the 
simulation runs, and on the initial value of A^*. Cur- 
rently, there are some observational constraint on the 
rate of mergers, which will be considered in a foUowup 
study. In this study, the results are merely normalized ar- 
bitrarily to match the final A^BH-A^gai relation observed 
today. 

The simulation "clock. " — The progress of time is not 
well defined in Monte-Carlo simulations, so it is useful 
to define merging cycles for the sake of keeping track of 
the simulation progress. Each full cycle is defined as be- 
ing complete after the number of merger events equals 
the number of galaxies present at the beginning of that 
particular cycle. Galaxies that are produced or merged 
retain their states for the following cycle. Because some 
galaxies may merge multiple times by being drawn re- 
peatedly, not all galaxies will be involved in mergers af- 
ter each full cycle. The exact definition of the simulation 
clock is unimportant, as it is only the relative number of 
major vs. minor mergers on average that determines the 
degree of convergence toward a linear A^BH-A^gai rela- 
tion, where a major merger is defined as having a mass 
ratio of at most 4:1. 

3.2. Replenishment Scenarios 

The replenishment scenario is one of the two simple 
ways considered to emulate the progress of galaxy evo- 
lution. Here, by definition, the rate of galaxy mergers 
equals the rate of galaxy number production. The way a 
galaxy is newly produced is by being selected randomly 
from an initial mass function, parameterized by a and 
A4* , which does not evolve with time. In contrast, galax- 
ies "grow" only by merging with another member in the 
galaxy pool existing at the time, and hierarchical merg- 
ing is the only avenue for mass growth. Therefore, as 
galaxies merge, the cumulative mass function does un- 
dergo evolution. However, the total number of galax- 
ies remains constant because of the 1:1 ratio of merg- 
ing : replenishment . 

Figures 3 and 4 illustrate the results for initial con- 
ditions P = (i.e. no A^BH-A^gai correlation) and 
(3 = 2 (steep A^BH--A^gai correlation), respectively, show- 
ing only a small subset of the data points. In each figure, 
two different initial mass functions, a = —0.5 (Figs. 3a, 
4a) and a = —1.5 (Figs. 36, 4&) are considered. In each 
of the Figures, the initial distribution of the A^BH--A^gai 
relation (or lack thereof) is shown with crosses. The open 
colored data points show the AlBH--A4gai development of 
galaxies that have undergone N„iaj > 5 major merger 
episodes, after 10 (blue triangles), 100 (green squares). 
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Fig. 3. — No initial correlation (/3 = 0) in the A^BH-A^gal relation, in the replenishment scenario. The black crosses represent the initial 
distribution of points, and the solid line shows the local A^BH-A^gal relationship from iHaring &: Rixl II2004I ') — it is not a fit to the data 
points. The colored data points represent objects that have undergone at least five major merger episodes after fO (blue triangles), 100 
{green squares), and 1000 (red circles) complete "merging cycles." The crosses are the primordial distribution, corresponding to the initial 
mass function. The shaded region illustrates the locus of all points after 1000 cycles; the density of points doubles with each contour level. 
The cumulative histograms after the corresponding merger sequences are shown below the data points, (a) An initial Schechter power law 
slope of a = —0.5. (fe) An initial Schechter powerlaw slope of a = —1.5. 




Fig. 4. — Similar to Figure 3, except for a steep initial correlation (/3 = 2) in the A^BH-A^gal relation. See Figure 3 for details. 
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Fig. 5. — Similar to Figure 3b, showing the effect of central-hmit 
tendencies with increasing number of major mergers for galaxies 
after 1000 merger cycles. The colored data points illustrate objects 
that have undergone 1 < N^^^j < 4, (blue triangles), 5 < N^^j < 
14 (green squares), A^maj > 14 (red circles) major mergers. The 
greyscale contours shows the locus of all the points. 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




log(M^,,/Mj 

Fig. 6. — Similar to Figure 3, except the BH mass is drawn from 
a Schechter law of a = —1.5 instead of a Gaussian distribution. 
See Figure 3 for details. 



and 1000 (red circles) merger cycles have transpired. 
These data points effectively illustrate the progress of 
the A^BH--A^gai evolution for objects that might be mor- 
phologically identified as early-type galaxies of each cy- 
cle. For clarity, the contour levels represent the locus of 
points after 1000 merger cycles, and the levels are spaced 
at multiples of 2 in density. The luminosity functions of 
the galaxy pool at the end of the merger cycles are shown 
in the lower half of each diagram in corresponding col- 
ors and locale in mass. Lastly, a linear reference line is 
overplott ed in the Figures w ith normalization given by 
7?,o=800 (jHaring fc Rix|[2004 ). and the simulations are 
scaled/shifted arbitrarily to match; it is not a fit to the 



data points. 

As shown in Figures 3 and 4, the convergence towards 
a tight linear relation is fairly quick. After five major 
merger episodes a linear relation starts to emerge regard- 
less of the initial conditions of the mass function or the 
form of the A^BH-A^gai correlation. One reason for this 
quick convergence is the central-limit behavior of major 
mergers which is shown in Figure 5, in which the increas- 
ing number of mergers is represented by different symbols 
and shades. The one notable case where the convergence 
toward linearity is slower than the other scenarios is Fig- 
ure 46, where the effect is only evident at 10^°-^ Mq or 
greater, even as the scatter has decreased markedly. In 
general, if the A^BH--A^gai correlation is steep initially, 
the tail at low mass remains steep after a large number 
of major mergers has occurred, even as the massive end 
converges toward linearity. 

Lastly, the qualitative convergence effects do not de- 
pend on the assumption about the distribution of BH 
mass at each galaxy mass. Figure 6 shows an example 
that is in direct analog to Figure 3b, except that the BH 
mass is instead drawn from a Schechter mass function 
with a = —1.5. 



3.3. Depletion Scenarios 

The other extreme of the merger simulations is to con- 
sider what effect galaxy depletion from a finite reservoir 
has on the A^BH-A^gai relation. Because the number 
density of galaxies builds up over time, the depiction 
scenario is expected to not be realistic. Nevertheless, it 
is useful for illustrating how the TWBH-A^gai convergence 
is affected by a different evolution in the mass function 
as compared with the replenishment scenario. 

The depletion scenarios are constructed by creating a 
large sample of 5 x 10^ objects, initially having no corre- 
lation between BH and galaxy masses (Figure 7) or with 
a /3 — 2 correlation between the two (Figure 8). The BH 
masses are assigned to the galaxies with a log-normal 
distribution of dispersion a = 2 centered around Equa- 
tion 2. In each scenario, galaxies are created to have 
initial mass functions of a = —0.5 (Figures 7a and 8a) or 
a = —1.5 (Figures 7b, 8b). Then, as galaxies merge, no 
new ones are created to replace them. As a consequence, 
the mass function evolves by growing in A^*, the num- 
ber density decreases, and a sharp truncation develops 
at low masses (see lower half of Figures 7 and 8). As 
the number of merging cycles increases, the scatter de- 
creases quickly and converges toward a linear relation, as 
illustrated by the solid line. Once again, as shown in Fig- 
ure 8 (especially 8b) , the convergence is much slower for 
steep a and steep (3 compared with other scenarios. And 
while the convergence trends are noticeable, because of 
a dearth of minor galaxies with which to merge at late 
times (red circles) , the slope is virtually "frozen in," and 
the subsequent convergence is due mostly to the central- 
limit theorem. 

4. DISCUSSION AND CONCLUSION 

This study has revisited the issue of how galaxy merg- 
ing may affect the AlBH--A^gai scaling relation from 
the standpoint of basic mass addition and statistics, 
thereby clearly isolating the merger cause from other de- 
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Fig. 7. — No initial correlation (/3 = 0) in the A^BH--Mgal relation, depletion scenario. The contours represent the initial distribution 
of points, and a solid line shows the local A^BH-A^gal relationship from iHaring fc Rixl (2004) - it is not a fit to the data points. The 
colored data points represent objects that have undergone at least 1 major merger episodes after 1 complete merging cycles {blue triangles), 
10 (green squares), and 14 (red circles). The cumulative histograms after the corresponding merger sequences are shown below the data 
points, a) An initial Schechter powerlaw slope of a = —0.5. 6) An initial Schechter powerlaw slope of o = —1.5. 




4 6 8 10 12 14 4 6 8 10 12 14 

log(Mgai/Mo) log(M^„/Mj 

Fig. 8. — Similar to Figure 7, except for a steep initial correlation (/3 = 2) in the A^BH-A4gai relation. See Figure 7 for details. 
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tailed physics that must otherwise affect galaxy evolu- 
tion. Through Monte-Carlo simulations, a tight, linear, 
■A^BH--A^gai correlation appears to emerge when galaxies 
have undergone five or more major mergers (along the 
entire tree, not just the main branch), and many minor 
ones, for practically all reasonable initial correlations be- 
tween A^BH and Algal, or a lack of one. The main reasons 
for these behaviors are seen to be the following: 

1 . The galaxy mass function decreases with increasing 
mass. 

2. Major mergers have a strong central-limit ten- 
dency, so that regardless of the initial A^BH--A^gai 
correlation, the scatter should decrease with an 
increasing number of events. While this ten- 
dency also acts on minor mergers, the drive toward 
smaller scatter is weaker because minor mergers oc- 
cur between galaxies that are vastly discrepant in 
both galaxy and BH mass as compared with major 
mergers, by definition. The corollary is that the 
steeper a correlation between A^bh (j/-axis) and 
Algal J the stronger the central- limit tendency for 
major mergers compared with minor. However, 
major mergers alone are not enough to cause the 
AleH-Algai relation to converge to linearity over 
time because the ratio A^Bn/A^gai is not changed 
much. 

3. Minor mergers are primarily responsible for causing 
the AIbh- Algal relation to converge toward a linear 
— that is, AIbh = LAigai — relation because the 
mass function of galaxies follows a Schechter pow- 
erlaw. Without minor mergers, the AlBH-Albuigo 
relation can be "frozen" to a slope that is not 
necessarily linear. This linear attractor causes a 
convergence toward a tighter AIbh- Algal relation; 
however, it is less efficient at reducing the scatter 
compared with the central-limit seeking tendency 
of major mergers, as shown in Figure 5. 

It is curious that galaxy merging itself might produce 
a linear AleH-Algai relation. However, a natural ques- 
tion that does arise is, "When is the merging statistics 
presented in this study relevant?" On the surface, it 
is easy to conclude that because the reasoning refers 
to a two component model it ought to apply to "dry" 
mergers, but perhaps not to a three component model 
involving stars, gas, and BH. Thus, the implication is 
also that it ought not apply to galaxies undergoing gas- 
rich mergers, that is, early cosmic history. However, it 
is not clear that such a skepticism is warranted. For 
example, in the entire discussion thus far, the abscissa. 
Algal, might just as well refer to Algai=A(stciiar+Algas, 
instead of just Alstoiiar- If BHs do not grow much by 
accretion and that the gas does not get removed from 
the definition of AJgai during mergers, then the A^bh- 
Aigai correlation can emerge from statistical merging. 
The argument holds true even if At gas transforms arbi- 
trarily into AJstoiiar, as long as the sum is conserved. In 
the limit where BHs do grow most of thei r mass dur- 
ing AGN accretion, as mig ht be implied bv lSoltanI (e.g. 
ll982D : lYu fc Tremaind fe.g. 120021 ) . so that ATWeHOcTUgas 
and AAIbh^AIbh, then the correlation between A^bh 



and Algal comes out by construction rather than by sta- 
tistical merging. However, statistical reasoning would 
still be a "supporting actor" to reduce the scatter and 
to forcibly steer the AlBH-A^gai relation in the preferred 
linear direction. Likewise, even if BH growth, or other 
physics (e. g. gravitational radiation, three body BH 
ejection - iMerritt et al.ll2004t IWlonteri fc Pemal 120051 : 
ICiotti et al.ll2007l and references therein), were a "heat- 
ing" source, that is, one that randomizes a tight linear 
A^BH-A^buigc relation, the linear and central limit attrac- 
tors would cause a re-convergence if galaxies continue to 
merge thereafter by both major and minor mergers. In 
summary, while it is entirely possible that the AIbh- 
A4 bulge relation has origins outside of basic statistics, 
galaxy merger statistics can still affect the final outcome 
of a AlBH-Algai correlation in both the scatter and the 
slope. In any event, statistical reasoning is a fundamen- 
tally robust explanation for why random galaxy merging 
does not corrupt a pre-existing AlBH-Albuigo relation, 
which is important to bear in mind in the context of the 
AlBH-Albuigo or A^BH-A^gai relation in a hierarchically 
forming universe. 

While the A^BH-A^buigo relation might have other 
origins, it is nonetheless interesting and revealing to 
follow through the consequences of statistical merg- 
ing. For instance, simple statistics naturally explains 
why black holes appear to correlate most strongly 
with galaxy bulges, rather than more generally with a 
galaxy as a whole, which mig ht include a stellar disk 
(jKormendv & Gebhardtl [2001 : bulge masses, assum- 
ing they were assembled through major mergers, have 
a stronger central-seeking tendency than disk galaxies, 
whose growth history might involve more minor merger 
events. As such, the A^BH-A^buigc relation is a special 
case of a more fundamental AlBH-A^gai relation. Revers- 
ing the argument, the observational fact that AIbh corre- 
lates most strongly with bulge masses, coupled with the 
central-limit theorem reasoning, implies that the merger 
trees of elliptical galaxies were more dominated by ma- 
jor merger events than were disk galaxies. Conversely, 
the fact that the scatter in the AlBH-Algai relation is 
observed to be much larger for disk dominated galaxies 
implies, statistically, that their progenitors, and progen- 
itors thereof, have undergone more minor mergers. 

The possibility that a more fundamental correlation is 
between Mbh and Algai (rather than Albuigc) also has 
practical implications for what slope and scatter would 
be measured by observations. First, because the slope 
changes with mass even for objects that experienced the 
same number of major mergers (e.g. Figures 3 and 7), 
the deviation from linearity and the intrinsic scatter will 
depend exactly on how the data are cut. Simply defining 
a sample of objects based on a mass selection cut will 
bias one's measurement of the slope and scatter. Fur- 
thermore, defining a sample based on morphology criteria 
may also implicitly preselect samples that have certain 
major vs. minor merger histories. Observationally, it is 
therefore crucial, when comparing intrinsic scatter and 
slope of the AlBH-Albuige relation to be specific about 
sample selection parameter space, morphology, bulge- 
to-disk ratios, or other criteria, lest the conclusions be 
caused by subtle but trivial selection biases. 

Another consequence of this thought experiment is 
that the ratio r=A^BH/A4buigc approaches an asymp- 
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totic value with time from having a smaller ratio in 
the past. O n the surface, this a ppears cont r ary t o 
the findings o f 'Peng et ^ (|2006allbl ): IWoo et all (|2006D : 
[Shields et all |2003, 200^ based on quasar host galaxy 
studies that the ratio F decreases over time. If the 
quasar host galaxy studies are correct and are not sig- 
nifican tly affected by biases pointed out by iLauer et al.1 
(l200l . then some other physics not considered here is 
responsible for causin g a decline in the normalization of 
r with time (e.g. see lCrotonl[200l iHopkins et al.l [20071 : 
iFontanot et al.l [20061 ). For instance, the abscissa is am- 
biguous about what mass A^gai corresponds. If gas mass 
is a significant fraction of a galaxy's mass, then form- 
ing stars out of the gas reservoir would decrease F over 
time, if the abscissa A^gai represents the galaxy's stellar 
bulge mass. Secular growth of galaxy bulges by accret- 
ing stars in galaxy disks would also decrease F, at the 
expense of increasing the scatter. Major mergers of pure 
stellar bulges, however, would not cause F to decrease 
over time. 

In hindsight, t he results of t his study co ul d have been 
anticip ated from [Islam et all (|2003l I2004D : ICiotti et all 
([20f)l . given that the initial conditions used in those 
studies are a special case of this one where the initial 
BH scatter ctbh (Figure 3a or 3b) (M. Volonteri 
and L. Ciotti 2 007, private communicati on). Just as rel- 
evant, [Croto3 (|2006( ): ICiotti et all (|2007[ ) show that once 
the A^BH-A^buige relation is in place, it is fairly imper- 
vious to being randomiz ed by galaxy merging. And be- 
cause [Ma^^i^n (|2003[ ) also uses realistic cosmological 
merger trees, they confirm that the arguments presented 
here ought to remain relevant. However, the reasons be- 
hind the AlBH--A4buigc convergence behavior are difficult 
to extract from previous studies because of the use of pri- 
ors, the use of identical BH seeds, the inclusion of other 
physics, and the focus on only the BH-bulge coevolution 
(i.e. major mergers). The latter, especially, is worth 
examining further, because the prior that one chooses 
about whether the BH correlates with just its bulge or 
with the entire galaxy can lead to differing interpreta- 
tions. 

In particular, one conclusion from [Islam et all (|2003l 
[2004[ ) is that the A^BH-A^buigo relation converges to a 
slightly non-linear slope oi /3 — 0.9; hence they reason 
that other physics, perhaps BH growth through accre- 
tion of gas, is required in order to increase the slope 
closer to linearity. The reasoning presented in the cur- 
rent study, however, would stipulate that linearity is an 
asymptotic outcome of mergers, but deviations from lin- 



earity come from the possibility that the low mass galax- 
ies have not yet achieved the asymptotic limit, because 
of a weaker convergence. At low masses, the slope de- 
viates from unity in either direction depending on the 
initial mass function of the galaxies (e.g. compare Fig- 
ure 3a with 3b) , on the mass cut of the study, and on the 
relative incidence of minor versus major mergers. 

An interesting consequence to consider is how the 
A^BH-A^gai relation might differ between high and low 
density galaxy clusters. However, one of the unrealis- 
tic side-effects of using Monte-Carlo simulations to de- 
termine merger rates is that the normalization of the 
A^BH--A^buigc relation is the same in all density regimes. 
This is because the normalization factor, F, depends only 
on the ratio of major to minor merger events which, in 
the Monte-Carlo universe, is not affected by a simple 
rescaling of the mass function. However, in the real uni- 
verse, the relative rates of major and minor mergers can 
change with density and, as such, may result in different 
normalization and scatter in the A^BH--A^gai relation. 

Lastly, because of the ambiguity in what Algai corre- 
sponds, depending on whether it refers to the total stellar 
mass, gas mass, dark matter halo mass, or a combination 
thereof, the degree of scatter and linearity would clearly 
differ, as a result of different initial mass functions and 
merger histories. Because the scenarios considered above 
depend on a linear addition of masses, the arguments 
therefore may not apply to gas masses that are not grav- 
itationally bound to a galaxy. These and other issues will 
be addressed in a future study, which will incorporate the 
use of realistic merger trees. 
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APPENDIX 

This appendix shows that the A^BH-A^buigc relation follows a central-limit-like behavior when galaxies undergo major 
mergers. Specifically, this means that if the initial parent distribution of progenitor Mbk that undergoes merging is, for 
simplicity, normally distributed about a mean BH mass fi, thus having a logarithmic dispersion cr(log(/iinit)), then a new 
distribution of AIbhs after merging will have a log-normal dispersion that scales as; cr (log (/imcrgc)) ~ cr(log(/iinit))/\/2 

First, the mean, /imcrgc, of the resulting BH distribution after two BHs merge from the initial parent distribution is: 

Mmcrgc (A^BH.l + A^BH,2) , (1) 

where A^bh.i and A^bh,2 are drawn from the same parent distribution for major galaxy mergers. Then, 



log (/Xmcrgc) = log {MbH,! + A4bH.2) , 



(2) 
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From propagation of errors the log-normal error is: a{\og{x)) = a{x)/x, then, 

cr(7WBH,l + A1bh,2) 



O- (log (/"merge)) = 



{MbH,! +XbH,2) ' 



(3) 



By definition, a distribution obtained by averaging the mass of merging BH pairs is a normal distribution with a mean 
of the initial parent distribution, /iinit: 



Minit 



Mmerge 



(4) 



Because Ad^i and A4^2 are drawn from a normalized distribution around a parent mean fj,, the new distribution of 

0-(/imerge) = CT {M,l + M ,2) is: 



a ill 



merged 



^(^init) 

V2 



X 2. 



Substituting Eqs. A4 and A5 into A3 yields: 



0-(l0g(Mmerge)) 



O-(Minit) 



Using the fact that, cr(log(/L()) = cr(/x)//x. Equation A6 becomes: 

a (l0g(/iinit)) 



<T(l0g(/iinerge)) 



V2 



(5) 



(6) 



(7) 



REFERENCES 



Adams, F. C, Graflf, D. S., Mbonyc, M., & Richstonc, D. O. 2003, 
ApJ, 591, 125 

Adams, F. C, GraflF, D. S., & Richstonc, D. O. 2001, ApJ, 551, 
L31 

Barth, A. J., Greene, J. E., & Ho, L. C. 2005, ApJ, 619, L151 
Barth, A. J., Ho, L. C, Rutledge, R. E., & Sargent, W. L. W. 2004, 
ApJ, 607, 90 

Borys, C, Smail, I., Chapman, S. C, Blain, A. W., Alexander, 

D. M., & Ivison, R. J. 2005, ApJ, 635, 853 
Boylan-Kolchin, M., Ma, C.-R, & Quataert, E. 2006, MNRAS, 369, 

1081 

Burkert, A., & Silk, J. 2001, ApJ, 554, L151 

Cattaneo, A., Blaizot, J., Devriendt, J., & Guiderdoni, B. 2005, 

MNRAS, 364, 407 
Cen, R. 2007, ApJ, 654, L37 

Ciotti, L., Lanzoni, B., & Volonteri, M. 2007, ApJ, 658, 65 
Ciotti, L., & van Albada, T. S. 2001, ApJ, 552, L13 
Croton, D. J. 2006, MNRAS, 369, 1808 

Croton, D. J., Springcl, V., White, S. D. M., Dc Lucia, G., Frcnk, 

C. S., Gao, L., Jenkins, A., Kauffmann, G., Navarro, J. F., & 
Yoshida, N. 2006, MNRAS, 365, 11 

Di Mattco, T., Colbcrg, J., Springel, V., Hcrnquist, L., & Sijacki, 

D. 2007, ArXiv c-prints, 705 

Di Mattco, T., Springcl, V., & Hernquist, L. 2005, Nature, 433, 
604 

Escala, A. 2006, ApJ, 648, L13 

Ferraresc, L., Cote, P., Dalla Bonta, E., Peng, E. W., Merritt, D., 
Jordan, A., Blakcslcc, J. P., Ha§egan, M., Mei, S., Piatek, S., 
Tonry, J. L., & West, M. J. 2006, ApJ, 644, L21 
Ferraresc, L., & Merritt, D. 2000, ApJ, 539, L9 
Filippenko, A. V., & Ho, L. C. 2003, ApJ, 588, L13 
Fontanot, F., Monaco, P., Cristiani, S., & Tozzi, P. 2006, MNRAS, 
373, 1173 

Gabasch, A., Hopp, U., Feulncr, G., Bender, R., Seitz, S., SagUa, 
R. P., Snigula, J., Drory, N., Appcnzellcr, I., Heidt, J., Mehlert, 
D., Noll, S., Bohm, A., Jager, K., & Ziegler, B. 2006, A&A, 448, 
101 

Gebhardt, K., Bender, R., Bower, G., Dressier, A., Faber, S. M., 
Filippenko, A. V., Green, R., Grillmair, C, Ho, L. C, Kormendy, 
J., Lauer, T. R., Magorrian, J., Pinkney, J., Richstone, D., & 
Tremaine, S. 2000, ApJ, 539, L13 

Giallongo, E., Salimbeni, S., Menci, N., Zamorani, G., Fontana, A., 
Dickinson, M., Cristiani, S., & Pozzetti, L. 2005, ApJ, 622, 116 

Graham, A. W., & Driver, S. P. 2007, ApJ, 655, 77 



Graham, A. W., Erwin, P., Caon, N., & Trujillo, I. 2001, ApJ, 563, 
Lll 

Granato, G. L., De Zotti, G., Silva, L., Bressan, A., & Danese, L. 

2004, ApJ, 600, 580 
Greene, J. E., & Ho, L. C. 2005, ApJ, 630, 122 
— . 2006, ApJ, 641, L21 

Haehnelt, M. G., & Kauffmann, G. 2000, MNRAS, 318, L35 

Haring, N., & Rix, H.-W. 2004, ApJ, 604, L89 

Ho, L. 1999, in ASSL Vol. 234: Observational Evidence for the 

Black Holes in the Universe, ed. S. K. Cheikrabarti, 157 — h 
Ho, L. C. 2007, arXiv:0707.3436 

Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., Robertson, 

B. , & Springel, V. 2006a, ApJS, 163, 1 

Hopkins, P. F., Hernquist, L., Cox, T. J., Robertson, B., & Krause, 

E. 2007, ArXiv Astrophysics e-prints 
Hopkins, P. F., Somcrville, R. S., Hernquist, L., Cox, T. J., 

Robertson, B., & Li, Y. 2006b, ApJ, 652, 864 
Islam, R. R., Taylor, J. E., & Silk, J. 2003, MNRAS, 340, 647 
— . 2004, MNRAS, 354, 427 

Kaspi, S., Maoz, D., Netzer, H., Peterson, B. M., Vestergaard, M., 

& Jannuzi, B. T. 2005, ApJ, 629, 61 
Kaspi, S., Smith. P. S., Netzer, H., Maoz, D., Jannuzi, B. T., & 

Giveon, U. 2000, ApJ, 533, 631 
Kauffmann, G.. & Haehnelt, M. 2000, MNRAS, 311, 576 
Kormendy, .1., & Gebhardt, K. 2001, in AIP Conf. Proc. 586: 20th 

Texas Symposium on rclativistic astrophysics, ed. J. C. Wheeler 

& H. Martel, 363-+ 
Kormendy, J., & Richstonc, D. 1995, ARA&A, 33, 581 
Lauer, T. R., Fabcr, S. M., Richstone, D., Gebhardt, K., Tremaine, 

S., Postman, M., Dressier, A., AUer, M. C, Filippenko, A. V., 

Green, R., Ho, L. C, Kormendy, J., Magorrian, J., & Pinkney, 

J. 2006, ApJ, 662, 808 
Lauer. T. R,., Tremaine, S., Richstone, D., & Faber, S. M. 2007, 

arXiv:0705.4103 

Li, Y., Hernquist, L., Robertson, B., Cox, T. J., Hopkins, P. F., 
Springel, V., Gao, L., Di Matteo, T., Zentner, A. R. , Jenkins, 
A., & Yoshida, N. 2006, ArXiv Astrophysics e-prints 

Magorrian, J., Tremaine, S., Richstone, D., Bender, R., Bower, G., 
Dressier, A., Fabcr, S. M., Gebhardt, K., Green, R., Grillmair, 

C. Kormendy, J., & Lauer, T. 1998, AJ, 115, 2285 
Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 
McLure, R. J., & Dunlop, J. S. 2001, MNRAS, 327, 199 
Merritt, D., Milosavljevic, M., Favata, M., Hughes, S. A., & Holz, 

D. E. 2004, ApJ, 607, L9 

Miralda-Escude, J., & KoUmeier, J. A. 2005, ApJ, 619, 30 



12 



PENG 



Nipoti, C, Londrillo, P., & Ciotti, L. 2003, MNRAS, 342, 501 
Onken, C. A., Fcrrarcsc, L., Merritt, D., Peterson, B. M., Pogge, 

R. W., Vcstcrgaard, M., & Wandel, A. 2004, ApJ, 615, 645 
Peng, C. Y., Impcy, C. D., Ho, L. C, Barton, E. J., & Rix, H.-W. 

2006a, ApJ, 640, 114 
Peng, C. Y., Impey, C. D., Rix, H.-W., Kochanek, C. S., Keeton, 

C. R., Falco, E. E., Lehar, J., & McLeod, B. A. 2006b, ApJ, 649, 

616 

Peterson, B. M., & Bentz, M. C. 2006, New Astronomy Review, 
50, 796 

Robertson, B., Cox, T. J., Hernquist, L., Pranx, M., Hopkins, P. P., 

Martini, P., & Springel, V. 2006a, ApJ, 641, 21 
Robertson, B., Hernquist, L., Cox, T. J., Di Matteo, T., Hopkins, 

P. P., Martini, P., & Springel, V. 2006b, ApJ, 641, 90 
Salviander, S., Shields, G. A., Gcbhardt, K., & Bonning, E. W. 

2006, New Astronomy Review, 50, 803 
— . 2007, ApJ, 662, 131 
Sclicchtcr, P. 1976, ApJ, 203, 297 

Shields, G. A., Gebhardt, K., Salviander, S., Wills, B. J., Xic, B., 
Brotherton, M. S., Yuan, J., & Dietrich, M. 2003, ApJ, 583, 124 



Shields, G. A., Salviander, S., & Bonning, E. W. 2006, New 

Astronomy Review, 50, 809 
Sijacki, D., Springel, V., Di Matteo, T., & Hernquist, L. 2007, 

ArXiv e-prints, 705 
Silk, J., & Rees, M. J. 1998, A&A, 331, LI 
Soltan, A. 1982, MNRAS, 200, 115 

Springel, V., Di Matteo, T., & Hernquist, L. 2005, MNRAS, 361, 
776 

Treu, T., Malkan, M. A., & Blandford, R. D. 2004, ApJ, 615, L97 
Treu, T., Woo, J.-H., Malkan, M. A., & Blandford, R. D. 2007, 

ArXiv e-prints, 706 
Vestergaard, M. 2002, ApJ, 571, 733 

Vestergaard, M., & Peterson, B. M. 2006, ApJ, 641, 689 
Volontcri, M., Haardt, F., & Madau, P. 2003, ApJ, 582, 559 
Volontcri, M., & Perna, R. 2005, MNRAS, 358, 913 
Woo, J.-H., Treu, T., Malkan, M. A., & Blandford, R. D. 2006, 

ApJ, 645, 900 
Yu, Q., & Tremaine, S. 2002, MNRAS, 335, 965 



